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To investigate the non- Fermi liquid behavior we consider the extended two-dimensional 
t — J model which includes additional hopping t". In the regime t, J <C t" we were able to 
solve the model analytically. It has a very rich phase diagram including antiferromagnetic 
(AF) insulator and AF strange metal with different kinds of pseudospin-singlet supercon- 
ducting pairings (p,d,g- waves). We also demonstrate a collective triplet excitation with 
energy below the superconducting gap. 



1 Introduction 

It is now widely accepted that superconductivity of cuprates is closely related to their 
unusual magnetic properties, and it is increasingly clear that magnetic pairing is the most 
realistic mechanism of cuprate superconductivity. However the mechanism of pairing 
as well as other unusual properties are far from completely understood. The problem 
has been attacked along several directions. First we have to mention the empirical or 
semi-empirical approach which allows one to relate different characteristics measured ex- 
perimentally. This approach is to a large extent based on the Hubbard model. For a 
review see article [p]]. In the low energy limit the Hubbard model can be reduced to the 
t — J model. Another approach to cuprates is based on numerical studies of the t — J 
model (see review 0). Our studies are also based on this model. We used the ordered 
Neel state at zero doping as a starting point to develop the spin- wave theory of pairing || . 
The method we used was not fully satisfactory, since it violated spin-rotational symme- 
try, nevertheless it allowed us to calculate from first principles all of the most important 
properties including the critical temperature, the spin-wave pseudogap and the low energy 
spin triplet excitations ||. 

A sharp collective mode with very low energy has been revealed in YBCO in spin 
polarized inelastic neutron scattering []5j, || [7]. A number of theoretical explanations 
have been suggested for this effect || |J, all of these are based on the idea that the 
system is close to AF instability. However all known explanations use some uncontrolled 
approximations and assumptions. 
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In the present work we investigate close to half filling regime for the 2D t — J model, 
where it can be solved analytically without any uncontrolled approximations. It can be 
done for the region of parameters where long-range AF order is preserved under doping. 
This is the regime where non Fermi liquid behavior can be studied in detail. We analyze 
the superconducting pairing in this regime and consider the spin triplet collective excita- 
tion. It is demonstrated that close to the point of AF instability energy of this excitation 
is very small. The excitation exists only at very small momenta. The idea of this work is 
somewhat similar to that of our previous paper || , however here we investigate different 
regime. 



2 Hamiltonian and single hole dispersion 

Let us consider a t — J — J" — V model defined by the Hamiltonian 

J- „ + n sr J 



H — t Cicr C jcr t C ia Cj 2Cr + 

(ij)a {ij 2 )a (ij) 



J ( SjSj — -riiUj j \ ;/,//, 



cJ CT is the creation operator of an electron with spin a (cr =f, [) at site i of the two- 
dimensional square lattice. The c\ a operators act in the Hilbert space with no double 
electron occupancy. The (ij) represents nearest neighbor sites, and (ij 2 ) represents next 
next nearest sites. The spin operator is Sj = ~ J2a,p H a a oiP c iPi an d the number density 
operator is rii = J2a c L c ?o-- m addition to the minimal t — J model (see Ref. 0) we have 
introduced additional next next nearest hopping t", and Coulomb repulsion V at nearest 
sites. Note that we do not introduce next nearest neighbor hopping t' (diagonal) because 
we do not need it for the purposes of this study. 

In the paper || we analyzed the model defined by the Hamiltonian ([!]) in the limit 
t, t" J. In the present work we consider limit 

t, J < t". (2) 

It is well known that the t — J model at half filling describes the Mott insulator. It is 
equivalent to the 2D Heisenberg model, and the ground state of the model has long range 
AF order. At small doping the holes are concentrated near the points (±7r/2, ±vt/2) where 
single hole dispersion has minima. In leading approximation the dispersion is of the form 
(we take energy at the minimum as a reference point) 

^ = /5(7^ + (7 k ") 2 ), (3) 
(3 w 0.8 x 8t" = 6At" 



i-'i 



7k = \{cosk x + cosfcj,), 7k = \{cosk x — cosk y ). Calculation of the dispersion (|3]) is 
straightforward because it is due to hopping within the same magnetic sublattice. Coeffi- 
cient 0.8 appears because of spin quantum fluctuations: 0.8 = 1 — 0.2, where 0.2 is the spin 
flip probability in the Heisenberg model. Alon with quasimomentum the hole in the AF 
background has an additional quantum number: pseudospin. We denote the hole creation 
operator by h^ a , where a = ±1/2 is pseudospin. The relation between pseudospin and 
usual spin is discussed in the paper JlO . 
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We will consider the case of very small doping, 5< 1, with respect to half filling (total 
filling is 1 — 5) . In this case all holes are concentrated in small pockets around the points 
k = (±7r/2, ±7r/2). Single hole dispersion (|3]) can be expanded near each of these points 

e k = (4) 

where p = (pi,p2) is deviation from the center of the pocket: p = k — ko, p\ is orthogonal 
to the face of the magnetic Brillouin zone, and P2 is parallel to the face (see Fig. 1). The 
Fermi energy and Fermi momentum for the holes equal ep ~ \n(35, pp ~ (ttS) 1 ^ 2 . 
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Figure 1: Magnetic Brilloun zone and single hole dispersion 



3 Hole-spin-wave interaction and instability of the 
Neel state 

Spin-wave excitations on an AF background are usual spin waves with dispersion u; q = 



2JJ1 — 7^ ~ V%Jq, at q « 1, see Ref. [1 1 for review. The hole-spin-wave interaction 



is well known (see, e.g. Ref. [ITS 



H, 



h.sw 



I^k.qUk+qj^kTttq + ^k+qT^/Jq + H.C.J, (5) 
k,q V J 

5f k ,q = 4t\/2(7 k D r q + 7k+qVq), 

where = Ck,-^ is the hole creation operator with pseudospin a, ct q and M are the spin 
wave creation operators for S z = =Fl, and £/ q = + \ and 

Vq = —sign(jq)^J-^- — | are parameters of the Bogoliubov transformation diagonalizing 



the spin- wave Hamiltonian, see Ref.[0]. Virtual spin wave emission gives a correction to 
the hole dispersion, see Fig. 2. However this correction is small 5e ~ t 2 /t" and therefore 
can be neglected compared with (H). 
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Figure 2: Spin-wave correction to the single hole dispersion. Solid line corresponds to the 
hole and dashed line corresponds to the spin wave. 

To describe renormalization of the spin wave under doping, it is convenient to introduce 
the set of Green's functions HTSfl 



£>«a(*,q) = -i(7K(*)at(0)]), (6) 

D a p{t,q) = -i(T[a q (*)/3_ q (0)]), 

D Pa (t,q) = -i(T[/3i q (t)a q (0)]), 

D^fcq) = -i(T[/3l q (0/3- q (0)]). 

In the present work we consider only the long-range dynamics: g ~ ^ ~ pj? < 1. In this 
limit all possible polarization operators coincide || P aa (u,q) = P a/3 (u!,q) = Pp a (u,q) = 
Ppp(u, q) = II(a>, q), where II (u;, q) is given by the diagram presented at Fig. 3. 




Figure 3: Spin-wave polarization operator. Solid line corresponds to the hole and dashed 
line corresponds to the spin wave. 

For stability of the system the condition (Stoner criterion) 

c^ + 2n(0,q)>0 (7) 



must be fulfilled ||14|| . Otherwise the Green's functions (|6|) would possess poles with 
imaginary uj. Considering holes as a "normal Fermi liquid" [TJ| one can easily calculate 
the polarization operator at q <C pp\ 11(0, q) ~ — At 2 \/2q / n (3 , Ref. [fl^| . Relatively weak 



pairing, which we consider below, does not influence this result. Then the condition of 
stability can be rewritten as 

= QAt" > ^. (8) 

71 J 

To provide stability of the AF order we have to choose 

t" > t" c w OAt 2 /J. (9) 

If t < J or t ~ J the stability condition is automatically fulfilled since in the present work 
we consider t, J <C t". However at t ^> J one can violate the condition (H). In this case 
we will assume that t" > t". If t" is close to t" c it is convenient to introduce the parameter 

v 2 = i-^ = (t"-t:)/t» (io) 



7T.jp 

as a measure of this closeness. The criterion (0) is proportional to this parameter. 
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4 Spin-singlet p-wave pairing caused by the short- 
range attraction 

It is not convenient to consider the superconducting pairing in the magnetic Brillouin zone 
with four half-pockets (see Fig. 1). Because of this we translate the picture to the shifted 
zone with two whole pockets, Fig. 4. We stress that this is question of convenience only, 
the representations are absolutely equivalent because of the translational invariance. 























-71 





Figure 4: Shifted zone with two whole pockets. Shadowed area corresponds to occupied 
hole states. 

There are two mechanisms for the superconducting pairing: short-range attraction 
and long-range attraction. First we consider the short-range effect. Attraction between 
holes at nearest sites (short-range) is due to the reduction in number of missing AF links. 
The value of this attraction immediately follows from eq. ([I]) 



U = J(SiSj - - + V) « -0.58 J + V. 



Strong enough Coulomb repulsion (V > 0.58J) kills this mechanism. In the momentum 
representation the interaction flTT| ) can be rewritten as 



H V —W ^2 7ki-k 3 ^k3t^k4|^k2|^kiT^ki+k2,k 3 +k4- 
ki,k2,k 3 ,k4 



:i2) 



For scattering inside a hole pocket the interaction is practically momentum independent 
because k x »s k 2 ~ k 3 « k 4 ps (n/2,n/2), and hence 7ki-k 3 ~ 1- Such interaction gives 
"s-wave pairing" with the gap without nodes at the Fermi surface. The value of the 
superconducting gap one can easily find using the results of papers (TO, O. This gives 



A = Ct"VSe^ 4U = Ct"V5e- 5t "/^ 8J - y \ 



(13) 



where C ~ 10 is some dimensionless constant. The solution is valid only if V < 0.58 J, for 
stronger Coulomb repulsion the pairing disappears. It is important to stress the peculiar 
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symmetry properties of the above pairing. This peculiarity comes from the presence 
of long-range AF order. As we already mentioned, the gap has no nodes at the Fermi 
surface and from this point of view it is "s-wave pairing". However we remind that we 
have considered the pairing in the shifted zone and in this zone it is not easy to classify 
the states by parity. For well defined parity we have to return to the magnetic Brillouin 
zone, so we have to translate the outside parts of the Fermi surface by the inverse vector 
of the magnetic lattice G = see Fig. 5. 




Figure 5: Translation from the shifted zone to the magnetic Brillouin zone. The su- 
perconducting gap has no nodes at the Fermi surface. The gap changes sign under this 
translation. 

The point is that under such translation the superconducting gap changes the sign 
as it is shown at Fig. 5. This property follows from the fact that the coefficient in the 
interaction QT2D changes sign under such translation: 7 kl _ k3+G = — 7ki-k 3 (for details see 
paper 

Thus in reality we have negative parity pairing which is usually called p-wave. The 
above consideration was relevant to the hole pocket centered at (tt/2, 7r/2). Similar con- 
struction is valid for another pocket centered at (vr/2, —tt/2). Existence of two solutions 
corresponds to the double degeneracy of the E-representation of the C± v group. Taking 
linear combinations of the single pocket solutions we find two degenerate solutions for 
the entire Brillouin zone with lines of nodes k x = or k y = well outside the Fermi 
surface. We would like to stress that we have considered the spin-singlet (more exactly 
pseudospin-singlet) pairing! This situation is very much different from the usual one when 
p-wave pairing implies spin triplet. We repeat that the peculiarity is due to the presence 
of long-range AF order. 

5 D- and g-wave pairings caused by the long-range 
attraction 

The long range attraction comes from the spin-wave exchange shown on Fig. 6. In this 
exchange the typical spin- wave momenta are q ~ pf ~ \/S, and hence the typical distances 
are r ~ 1/q ~ l/v / 5> 1. 
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Figure 6: Spin-wave exchange mechanism of attraction. Solid line corresponds to the hole 
and dashed line corresponds to the spin wave. The arrow shows the hole pseudospin. 

Similarly to the previous section, it convenient to consider first the pairing inside a 
hole pocket, say centered at (7r/2,7r/2), see Fig. 4. This pairing has been considered 
in detail in our previous work ||. It has been shown that for the case of "isotropic" 
dispersion (|J) the only solution is the one with a single node line in the pocket. The gap 
at the Fermi surface (€f = \n[35) is of the form 

A(0) = A o sin0, (14) 
A = Ce F e~^ 2t2 « 10Ct"6 e ~ WJt " ' l \ 

where sin0 = P2/PF, v\ = Pi + P21 an d C ~ 1 is some constant. 

The eqs.(|i~4"D describe pairing within a single pocket of the shifted zone. Translation 
of this solution to the magnetic Brillouin zone is shown at Fig. 7. This is absolutely 
identical to what we did in the previous section (change of sign at the translation). 




Figure 7: Translation from the shifted zone to the magnetic Brillouin zone. The super- 
conducting gap has line of nodes. The gap changes sign under the translation. 

There are effectively two pockets in the Brillouin zone, see Fig. 4. Taking symmetric 
and antisymmetric combinations between the pockets, we get the d- and g-wave pairings 
respectively. The symmetries of the corresponding superconducting gaps are shown at 
Fig. 8. It is clear that the d-wave belongs to the Bi representation of the C± v group and 
the g-wave belongs to the A 2 representation. 

Both solutions originate from ([14]), therefore they are close in energy. Nevertheless 
the constant C in eq. ([TJ|) is smaller for the g-wave. This is the price for additional 
lines of nodes (k x = and k y = 0). The above consideration did not include short range 
interaction (|l2"l) . This is absolutely correct for g-wave pairing which is not sensitive to the 
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Figure 8: Symmetry of the superconducting gap corresponding to the d- and g-wave pair- 
ings. 

interaction (|T2|) at all. However the d-wave is sensitive. Therefore at V < 0.58 J the d-wave 
pairing is enhanced because of fli"2|), while, on the contrary, at larger Coulomb repulsion 
V > 0.58J the d-wave is suppressed and can even disappear. To avoid misunderstanding 
we stress that in the limit under consideration (t" ^> t, J) the short range interaction 



(p^) is too weak (even at V = 0) to produce d-wave pairing without spin-wave exchange. 
However the short-range interaction influences the dimensionless constant C (see eq. (|14|)) 
which arises in spin-wave exchange mechanism. 



6 The phase diagram 

The phase diagram of the model under consideration is given on Fig. 9. To be specific we 
present the case of the not too strong Coulomb repulsion at the nearest sites: V < 0.58 J. 
At stronger V the p-wave superconductivity disappears, see eq. (|13|). Comparing eqs. 
(|13|) and (|14D we see that the p-wave pairing is stronger at t < t c , while at t > t c the 
d-g-wave pairing dominates. At V = the critical value is t c ~ J. In the p-wave phase 
the gap, as well as the critical temperature, is proportinal to square root of the hole 
concentration: A ~ T c oc \^6. But in the d-g-wave phase they are proportional to the 
first power of concentration: A ~ T c oc 5 

According to eq. (0) at t < t C N ~ 1.6\/t"J, the long range AF order at zero tempera- 
ture is preserved under doping, so we have coexistence of the superconductivity and the 
Neel order. At t > t c ^ the Neel order is destroyed by the doping and one gets a transition 
into the quantum disordered phase. However as soon as the magnetic correlation length is 
larger than the superconducting correlation length the mechanism of pairing is valid and 
one still has the d-g-wave superconductor. At a temperature higher than the critical one 
the system behaves as a metal with very strong scattering of mobile holes on spin-wave 
excitations. Following the tradition we call this state "strange metal". 
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Figure 9: The phase diagram of the extended 2D t — t" — J — V model. t c is transition 
point from the p-wave to the d-g-wave superconductor. t cN is transition point from the 
Neel state to the spin quantum disordered state. 



7 The spin-wave collective excitation 

We will see that the spin-wave collective excitation has nontrivial behaviour in the vicinity 
of the quantum phase transition from the Neel to the disordered phase. Therefore we 
study this excitation only in the d-g-wave superconducting phase at T = 0. The energy 
spectrum and Bogoliubov parameters are given by the usual BCS formulas 



u k ,v k 



(15) 



2 V E k 



with gap Ak from eq. The spin wave polarization operator due to mobile holes 

is given by diagram on Fig. 3 plus a similar diagram with anomalous fermionic Green's 
functions. Straightforward calculation gives (see e.g. Ref.H]) 



2(E k + £ k _| 



k,kn 



q ^ 2 - (E v 



E: 



k+qj 



' 2 2 
M k f k+q 



WkfkWk+q^k+q 



(16) 



This equation includes summation over pockets k = (tt/2, ±7r/2). In these pockets the 
vertex (|) is (?k ,q ~ 2 5 / 4 t(g :E ± q y )/ ' ^/q. Let us consider the case of very small momenta 
and frequencies: v F q < A , and uj < A . In this limit one can put q = in eq. (|T6| ) 
everywhere except at the vertex and therefore the polarization operator can be evaluated 
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n ^ q) =-w( 1+ W 



analytically 



Note that the the imaginary part is nonzero even at uo < 2Ao because the gap fll4|) has a 
line of nodes. Any of the Green's functions (§) have a denominator cu 2 — a; 2 — 2co> q n(u;, q), 
see e.g. Refs. [|T3| , |4]]. The zero of this denominator gives the energy and width of the 
spin-triplet collective excitation. Using eqs.(|TTD and (|T0) we find 



Oq = TfjJq, (18) 
7T 1 - 7] 2 UJ^ 

T " = 8 — A^' 

In essence this is the renormalized spin-wave. Far from the point of AF instability the 
parameter 77 ~ 1 , therefore the renormalization is relatively weak and the decay width 
is small. The situation is different when approaching the point of instability t — > t c ^ « 
l.§Jt" / J. Here, according to eq. (|TDD , 77 — > and therefore the energy of the renormalized 
spin wave is much smaller than the energy of the bare spin- wave, o^/ujq = 77 <C 1. 
Moreover this collective excitation exists as a narrow peak only at very small q, when 

7Tuj8r)A < 1. (19) 

At higher q the width is larger than its frequency because of decay to particle-hole excita- 
tions. We stress that the closer to the point of instability, the smaller is 77, and therefore 
the smaller is the region of q where the excitation exists. 



8 Conclusions 

We have considered a close to half filling t — t" — J — V model at t" ^> t, J. We restrict 
our consideration to the case of small doping 5 1. It is demonstrated that at t < 
t C N ~ l.6Vt"J the Neel order is preserved under the doping, and at t > t cN the order is 
destroyed and the system undergoes a transition to the quantum spin disordered phase, 
see phase diagram at Fig. 9. 

If the hole-hole Coulomb repulsion at nearest sites is not too strong (V < 0.58 J), then 
at small t the model has psedospin-singlet p-wave superconductivity. As t increases, at 
the some point t c (at V = the critical point is t c ~ J) the system undergoes a phase 
transition from the p-wave to the d-g-wave superconductor, see Fig. 9. Which state is 
realized (d- or g-wave) crucially depends on the Coulomb repulsion V. If V is small the 
d-wave is preferable, while at larger V the g-wave superconductivity is realized. 

In the Neel state we found the collective spin triplet excitation (renormalized spin 
wave). In the vicinity of the quantum phase transition to the spin disordered state the 
excitation exists as a narrow mode only at very small momenta and its energy is substan- 
tially below the energy of the bare spin wave. 



10 



References 

[1] A. V. Chubukov and D. K. Morr, Phys. Rep. 288, 355 (1997) and reference therein. 
[2] E. Dagotto, Rev. Mod. Phys. 66, 763 (1994). 

[3] V. V. Flambaum, M. Yu. Kuchiev, and 0. P. Sushkov, Physics C 227, 267 (1994); 
V. I. Belinicher, A. L. Chernyshev, A. V. Dotsenko, and O. P. Sushkov, Phys. Rev. 
B 51, 6076 (1995). 

[4] O. P. Sushkov, Phys. Rev. B 54, 9988 (1996). 

[5] J. Rossat-Mignot et al, Physica C 185-189, 86 (1991). 

[6] P. Dai et al, Phys. Rev. Lett., 70, 3490 (1993); 77, 5425 (1996). 

[7] H. F. Fong et al, Phys. Rev. Lett., 75, 316 (1996); 78, 713 (1997). 

[8] E. Dernier ans S. C. Zhang, Phys. Rev. Lett., 75, 4126 (1995); D. Z. Liu et al, ibid. 
75, 4130 (1995); I. I. Mazin and V. M. Yakovenko, ibid. 75, 4134 (1995); V. Barzykin 
and D. Pines, Phys. Rev. B 52, 13585 (1995); G. Blumberg et al, ibid. 52, R15741 
(1995); N. Bulut and D. J. Scalapino, ibid. 53, 5149 (1996). 

[9] D. Marinaro and O. P. Sushkov. Phys. Rev. B 58, 14934 (1998). 

[10] O. P. Sushkov, G. A. Sawatzky, R. Eder, and H. Eskes, Phys. Rev. B 56, 11769 
(1997). 

[11] E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 

[12] C. L. Kane, P. A. Lee, and N. Read, Phys. Rev. B 39, 6880 (1989). 

[13] J. Igarashi and P. Fulde, Phys. Rev. B 45, 12357 (1992). 

[14] O. P. Sushkov and V. V. Flambaum, Physica C 206, 269 (1993). 

[15] "Normal Fermi liquid" is a somewhat misleading term because we consider a system 
which is not a normal Fermi liquid from the common point of view. In this case "Nor- 
mal Fermim liquid" means that we consider an ideal gas of Fermions with dispersion 
(0). At zero temperature all states inside the semicircles at the Fig. 1 are occupied. 

[16] M. Randeria, J.-M. Duan, and L.-Y.Shieh, Phys. Rev. Lett., 62, 981 (1989). 

[17] M. Yu. Kuchiev and O. P. Sushkov, Phys. Rev. B 53, 443 (1996). 

[18] M. Yu. Kuchiev and O. P. Sushkov, Physica C 218, 197 (1993). 



11 



